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Abstract 

We analyze the energy dependence for two types of parity-non-conserving (PNC) asymmetries 
in the reaction 'jD np in the near-threshold region. The first one is the asymmetry in reaction 
with circularly polarized photon beam and unpolarized deuteron target. The second one corre- 
sponds to those with an unpolarized photon beam and polarized target. We find that the two 
asymmetries have quite different energy dependence, and their shapes are sensitive to the PNC- 
meson exchange coupling constants. The predictions for the future possible experiments to provide 
definite constraints for the PNC-coupling constants are discussed. 
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I. INTRODUCTION 



For more than forty years, the parity non-conservation (PNC) in nuclear processes at- 
tracts attention as a unique tool for studying the strangeness conserving {AS = 0) weak 
nucleon-nucleon interaction defined by nontrivial interplay of the weak quark-quark inter- 
action and the QCD-dynamics of composite hadrons at short distances Q^Ql- Most of the 
present theoretical studies of the parity non-conservation in nuclear processes are based on 
the finite-ran|;e vr, u and p-meson exchange potential of Desplanques, Donoghue, and Hol- 
stein (DDH) [3]. Using the symmetry consideration and the constituent quark model, DDH 
found the "reasonable range" and the "best values" of the the PNC meson-nucleon coupling 
constants. Their predictions are related to the theory of the weak interaction. Thus, the 
"best values" of the tcNN coupling for the Cabibbo and Weinberg-Salam models correspond 
to /i^ ~ 0.2 and 4.6 (in units of 10"'^), respectively. The predictions for the vector meson- 
nucleon weak coupling constants are also "theory-dependent", but this dependence is not so 
strong. In case of the charge-current theory, the transition m — > s responsible for the nNN 
interaction is suppressed by tan^ Oq — 0.05 {6c is the Cabibbo angle) as compared with the 
other transitions. This results in strong reduction of /i^. The neutral-current theory is free 
from this suppression which leads to a large value of h^^. The value of depends also on the 
non-perturbative QCD-dynamics of interacting mesons and baryons. The predictions based 
on the Skyrmion model I^J, the QCD-sum rule ||5|, the soft-pion approximation 6], and the 
quark model with the A degrees of freedom [3] give the value of /ivr = 0.8 ~ 3 x 10^'' which 



br the Weinberg-Salam model), being 
for the review of these estimations ) . 



is in the "reasonable range" of the DDH prediction ( 
smaller than the corresponding "best value" (see Ref. 

Analysis of the available data from the nuclear PNC-experiments suggests that the 
isoscalar PNC nuclear forces dominated by the p and (X'-meson exchange are comparable 
with the DDH "best values", whereas the isovector interaction dominated by the vr-meson 
exchange is weaker by a factor of 3 21. For example, the measurement of the circular polar- 
i.«o„ of .he p«o. e...ed .o. F ..1. . t.e < < 1.8 (xlO^, bi- 

However, this constraint is in disagreement with the recent analysis of the ^^'^Cs anapole 
moment performed in Refs. [8|, . Quite different theoretical approaches result in sim- 
ilar conclusions: for adequate description of the data on the anapole moment, one needs to 
use which is about a factor of 2 greater than the DDH "best value" h^^'^^ ~ 4.6 ■ 10"''. 
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These experimental situations mentioned above impel the new measurement and theoretical 
studies to resolve subsisting inconsistencies. 

The studies of the PNC-transitions in the nucleon-nucleon are very attractive because 
the two-nucleon wave functions are known reasonably well. The reactions ^ np are 
particularly important. Up to now, great efforts have been devoted to analyzing the thermal 
neutron capture by proton in the reactions with unpolarized and polarized neutron. In 
this first case, the circular polarization of emitted 2.23-MeV photons is analyzed. The 
experimental value \PJ = (18±18)x 10~^ Il3|| is consistent with the theoretical estimations 

- n M n 

\P^\ = (1.8 ~ 5.6) X 10 ^ jlJ, ll^ lla]. But, poor accuracy does not allow to obtain 
any definite conclusion about the strength of the PNC-forces. In the second case, the 
subject of study is the spatial asymmetry of emitted photons. The experimental value 
of = (6 ± 21) X 10~^ |17] is again too crude to check the theoretical predictions of 
~ 5 X 10~*, ( see i.g. Ref. [isj for reference and quotations). At present, a new 
PNC-asymmetry measurement for the radiative neutron-proton capture is in preparation 
at LANSCE pjf in order to reduce the experimental error of A^. 

Different aspects of parity non-conservation in deuteron electro-disintegration were ana- 



lyzed in Refs. 



3, 
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22| |. However, the nuclear PNC-efFect in this reaction is found to be 



insignificant compared to the contribution of the 7 — Z-boson interference of the individual 
nucleons j^]. 

With the advent of the high-intensity polarized photon beams, investigation of PNC- 
effects in the •yD — > np reaction becomes very important 2^. It is clear that in this case one 
can obtain complementary information on the PNC-interaction. Moreover, one can study 
the dependence of the PNC-asymmetries as a function of photon energy (contrary to the 
radiative np-capture, where the photon energy is fixed: ^ 2.23 MeV). This allows to get 
additional information which might help to reduce the ambiguity induced by uncertainties 
of the parity-conserving NN-forces at short distances. Thus, for example, the constraints on 
the PNC meson exchange coupling constants are usually obtained from compilation of the 
data extracted from the different experiments P|. This analysis includes a model-dependent 
estimation of the PNC-matrix elements in quite different objects like two- and few-body 
systems, light and heavy nuclei with their own assumptions and approximations. The energy 
dependent asymmetries in the 'jD —>■ pn reaction allow to give the similar constraints using 
only one simplest nuclear system. 
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In this paper, we discuss two PNC-asymmetries. One is the asymmetry Arl in deuteron 
disintegration in the reaction with circularly polarized photons and unpolarized deuteron. 
This asymmetry is mainly defined by the AI = 0,2 PNC-interaction and is equal to 
at ^ -Ethr; where i^thr is the threshold energy. The second one is the deuteron spin 
asymmetry Ad in reaction with unpolarized beam and polarized deuteron target (polar- 
ized along-opposite to the beam direction). It depends also on the isovector A J = 1 PNC- 
interaction, and therefore may be used for examining h^r. The A/j^-asymmetry was analyzed 



previously in Refs. [2J, |25|, |26|. In Refs. j2J, |25|, the calculation has been done only with 



repulsive hard-core NN-potentials which seems to be obsolete compared to the more sophis- 
ticated realistic potentials with soft repulsion. Energy dependence of Aj^l in the region 
- Ethr ~ 0.5 - 5 MeV was skipped. In Ref. |2J], the contribution of the PNC-vrA^A^- 
transition was completely ignored. On the other hand, it was included in Ref. [2^, and the 
extraordinarily big contribution of the weak ttNN transition to Aju at E^— E'thr = 1 ~ 30 
MeV has been reported. This result was used by other authors (cf. e.g. jaj27|) to discuss 
a possibility for extracting /i^ from the Aju-asjmmetrj. However, in Ref. j2a], it is shown 
that the consistent description of all transitions defined by the spin-conserving AI = 1 inter- 
action results in their mutual cancellation which is a disadvantage of using A^l as a tool for 
studying the weak ttNN transition. In Ref. 26[, the PNC-asymmetry is calculated on the 
basis of zero-range approximation where the short-range behaviour of the proton-neutron 
wave functions is modified phenomenologically, and therefore this result may be considered 
as a raw qualitative estimation. The PNC-asymmetry, An, is analyzed in Ref. |3| within 



the same model as given in Ref. [2q and therefore its result remains at very qualitative level. 
In our study, we use two realistic NN-potentials. One is the Paris potential 29, 3^ with 



soft repulsion at short distances and another is the Hamada- Johnston (HJ) potential 
with hard core repulsion. The long-range meson-exchange part of the NN-interaction in 
these potentials coincides, and the difference appears at short distances. Our results with 
the Paris potential may be useful as a prediction for future possible experiments, because the 
Paris potential was designed specially for proper description of the short range phenomena. 
The results with the HJ-potential are rather illustrative, and we show them in order to link 
our calculation with the previous works and to show explicitly the effect of the short-range 
correlation as an example of the extreme hard repulsion. 

In calculations of the PNC-asymmetries, the usage of models motivated by QCD (i.g., the 
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effective chiral perturbation model (ChPT) 32,|33|) seems to be interesting and important. 



However, the present status of ChPT allows to use it on 



y for the processes dominated by the 



long-range A/ = 1 PNC forces (like A^-asymmetry j32]) which is not enough to be applied 
for our case where the short-range A/ = 0, 2 transitions are important. Therefore, we restrict 
to perform the present calculation only in the framework of the potential description. 

This paper is organized as follows. In Section II, we define observables for the regular 
Ml and i?l-transitions. The formula for the PNC-interactions and expressions for the odd- 
parity admixtures are given in Section III. In Section IV, we discuss the results and report 
some predictions for the future experiments. The summary is given in Section V. 

II. REGULAR TRANSITIONS 

Near the threshold with < 10 MeV, where is the photon energy, the deuteron 
disintegration np is dominated by the Ml transition D ^Sq and the £'l-transition 

D ^Pj- The amplitudes of these Ml and El transitions read: 

Ta(M1) = j dr^MS + + lp)[n x ex\^i, (1) 

Tx{El) = J drrfve^i;,, (2) 

where k = nfc is the photon momentum, £\ is the photon polarization vector, A is the 
photon helicity, M is the nucleon mass, /i^ = /^p + /Xn = 0.88 and /ii, = /Up — /x„ = 4.71 are 
the isoscalar and isovector nucleon magnetic moments, respectively; e is the electric charge 
a = e'^/AiT = 1/137; r is the proton-neutron relative coordinate: r = Tp — r„, Ip is the proton 
orbital momentum: Ip = —iTp x Vp = —ir x V/2 = 1/2; tpi and ipf are the proton- neutron 
wave functions in the initial and final states, defined in the obvious standard notations as 

ui{r) 



ij, = J2{^fila\lM,)Yi,{r)x, 



Mi ^ ' 

(2S+1 ; 



i^f = 4vr^(//i.a|JM^)y,;(p)ll,(r)x. /^ ^ (3) 

lsfJ,(T 

where uo{r) = u{r) and U2{r) = w{r) are the radial deuteron s- and d- waves, respectively, 
and u{^^~^^Kj : pr) {K = S,P ...) is the radial continuum wave function. The spin operators 
S and S in Eq. ((H) are defined as 

S = ^(cTp + '^n), S = ^(o-p - cr„). (4) 



The upper and lower signs in Eqs. (HJ and © correspond to the photon absorption or 
emission, respectively (iq . 

In the following consideration, the regular and PNC - transitions from the np - bound 
(with the radial wave function u^) to the np ^Pj scattering states (with the corresponding 
radial wave function uj) will appear. Our analysis shows that these radial integrals at 
considered energies are not sensitive to J, therefore we can use " degenerated" approximation, 
where uj is calculated with the central forces. The reason of week sensitivity of the radial 
integrals to J is that the dominant contribution to the radial integrals comes from relatively 
large distances, where mo,i,2 are close to each other because the phase shifts for different 
states at < 10 MeV are rather small: \6j\ < 4 degrees. Small distances with r < 0.5 fm, 
where uj are really different, does not contribute in the integral because uj are small, and 
because of strong suppression from uoir) (or ruoir)). Direct numerical calculation shows 
that for < 10 MeV the validity of this approximation is better than 4-5% which is quite 
reasonable. This approximation allows to express the corresponding matrix elements in a 
very transparent form useful for qualitative analysis. But this approximation can not be 
used for calculation of the odd parity admixtures. In this case, the spin-orbital and tensor 
parts of NN potentials have to be taken properly into account. 

The regular Ml and El transition amplitudes expressed through the radial proton- 
neutron wave functions have the following form 

T,(M1) = -AiV^/0, 5_,M„ I'm = J u*CSo : pr) u{r) dr, (5) 
T^{E1) = iN^m Yl (^f'^^\JMf)Y*^ip) [(^^A^.Af/^- v^(2mla|lM,)(2mlA|l/i)/|](p) 

/»^/..eP.:p,.),.(,.),.*.. Il^j u-i'P.r.rMr)r,r, (7) 

where u{r) and w{r) are the radial deuteron s- and d- waves, respectively, and p is the 
proton momentum in c.m.s. 

The normalization factor in Eqs. (0) and (jUj) reads 

= 8 
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The total cross section is related to the amplitudes Tx as 



a 

1 



$^(|T.(M1)|2 + \T,{El)\^), (9) 



ITaP = —j d%\Tx\\ (10) 
where Mj is the deuteron spin projection and 

^E^M^ = i/^r+^i/ir. (11) 

AA^, 

In the following, we will assume the average of Eq. (jlOj) in all quadratic forms of TaT^ which 
define the observables for the case when the angular distribution of the final nucleon is not 
fixed and skip the symbol "overline", for simplicity. 

The wave functions for the deuteron bound state and the np-scattering states are calcu- 
lated using the realistic nucleon-nucleon potentials for two extreme cases: potential with soft 
short-range repulsive core (Paris potential 0, 0|) ^^"^ potential with hard-core repulsion 
(Hamada- Johnston (HJ) potential 

Figure Q shows the result of our calculation for the total cross section of the — >■ np 
reaction as a function of the energy excess: Ai?^ = — -Ethr? where i^thr is the threshold 
energy Ethr = e(l + e/2(Mp + M„ ~ ^ — ^i and e = 2.23 MeV is the deuteron binding 
energy, together with available data [sj, S|- "T^^ result for the Paris potential is shown 
in Fig. ^(a), where each contribution from Ml and El transitions is also displayed. The 
difference between the Paris and HJ-potentials in the total cross section does not exceed 5% 
and disappears at ^E^ —>■ (cf. Fig. ^(b)) because the main contribution into the radial 
integrals of Eqs. and (0) at small AE^ comes from the relatively large distances with 
r ^ 1 fm, where the rap-wave functions calculated for all the realistic potentials are close 
to each other. This result is in agreement with those of the previous calculations performed 
with various realistic potentials (cf. Ref. 0| for references and quotations). 
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III. PNC-INTERACTION AND PARITY ODD ADMIXTURES 



The short range PNC potential is expressed in terms of p, uj and vr exchanges and has the 



following form jS, 

^PNC = 



37| 



M 



X (S{V,/p(r)} + (l + Xp)OV/,(r)) 



Ti X T2 



SV/,(r) 



where 



+ 



M 



:S{V,/4r)} + (l + x.)^^V/4r)) 



+ 2^i(^i^-^2)S{V,/.(r)} 



y2M 1 



Ti X T2 



(12) 



/^(r) ~ /p(r) 



47rr ' 



47rr 



O = -[cTi X 0-2]. 



(13) 



For the strong nucleon- meson coupling constants Qi and Xii ^6 use commonly accepted 



values 0: Qp = 2.79, = 8.37, g^, = 13.45, Xp = 3.71 Xuj = 
nucleon coupling constants, hi, are taken as the "best value" of Ref. for the Weinberg- 
Salam model. The sensitivity of the observables to h^^ will be discussed separately. For 
convenience. Table 1 shows all parameters used in the present work. Parity-odd admixture 
states ip to the deuteron wave functions and np-scattering states are defined in the first 
order of perturbation theory in terms of Schrodinger equation 



-0.12. The PNC meson- 



[E - Hpcli^ = VpNcV', 



(14) 



where Hpc is the parity-conserving Hamiltonian and VpNc is the parity-violating two-body 
potential. For the odd-parity ^Pi admixture in a deuteron with J = 0, we have the following 
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expression: 

uCPi : r) 



« ^iM,(r)Xoo. 

—7^ I dr 91 (-e;r,r ) 



t=uj,p 



V3 
-V2 



u{r-') 



-3h' 



-X./;(r')+2/.(r')(|; + ^) 



w{r') 



(15) 



where XsSz nucleon spin function, gj^^E; r, r') is the Green function of the radial 

Schodinger equation for the np-system with the orbital momentum / = 1, isospin / = 0, spin 
S* = and the energy E = — e. 

The odd-parity ^Pi admixture with / = 1 is dominated by the 7r-meson exchange weak 
interaction. Nevertheless, for completeness we also include the contributions of the p and 
ci;-meson exchanges for AJ = 1 transition. The net expression for the ^Pi admixture reads 



ui^Pi : r) 



;^(l/il(j|lMi) 



ui^Pi : r) 



{1(7 



= ^y"rfr'<(-e;r,r') I /:(r')-v^^7p/.';,/;(r'^ 
- V2{gX^ - g,hl) [(/;(r') + 2/p(r')(^ - ^))n(rO 
+;^(/;(0+2/p(r')(|; + ^))^(r' 



'u(r') H — —w{r') 
v2 



(16) 



Figure El (a) shows the odd-parity ^Pi and ^Pi admixture in the deuteron wave function 
for the Paris (solid curves) and HJ (dashed curves) potentials. The main difference between 
the two potentials appears at short distances. In case of the HJ potential, all wave functions 
vanish in the core- region with r < Tcore (^core = 0.48 fm). This results in a sizeable suppres- 
sion of ^Pi-admixture because the "form factors" /^(r) in Eq. pSj) decrease sharply with r. 
The function hTr{r) decreases more slowly. Therefore, the ^Pi-admixture is not so sensitive 
to the choice of the potential model. 

Analysis of the odd parity component in the continuum np-states shows that at E^ < 
10 MeV, the dominant contribution to the considered asymmetries comes from the ^Pq 
admixture to the ^Sq state, from the ^Sq admixture to the ^Pq state, and from the ^Si and 
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^Di-components of the ^Pi-state. They are defined as follows 



An u{ Pq : pr) 



pr 



(17) 



uCPo-.pr) =- J2 ^9^i j dr'g\\E-r,r') 



t = p,UJ 



X 



(2 + X.)/;(rO + 2/.(r')(|;-^) 



uQSq : pr'), 



^('5o) 
u{}Sq : pr) 



XooZ^^lmlP) 




pr 



(18) 



I dr'gl%E;r,r') 



ui^Si : pr) 



V3 



u(^Si : pr) 
pr 



Xrar')-2Mr'){^ + ^] 



u{^Po : pr') 



XlM 



(19) 



dr'g',\E;r,r') Ig^Kf^r') - V2g^h'\fy) 



+ ^2{g^hl - g,hl) (/;(/) + 2/,(r') ( A + 1) 



m(^Pi : pr') 



ui^Di : pr) 



^^^ ^(^^1 :H ^(2;,la|lM^)y,,(r)xi.^y;^(p), 



pr 



(20) 



/iff 



rfr'^7°'(E;r,r') g^K f^') - V2g,h'lfy) 



d 2 



+ V~2{g^hl - g,hl) (^/;(r') + 2/,(r')(— - -) 



u{^Pi : pr'), 



3"- 



where E = p'^/M. The Green functions g{E; r, r') in Eqs. psp - psp are expressed through the 
regular and irregular solutions of the corresponding Schrodinger equations in the standard 
way. For the ^Si and ^Di states we use their spectral representation 



MgnE;ry) = ^J^^ + ^ I dk^'""^ ^ ^^^^^^^^ ^ '^'^ 



(21) 



E-E^ 

with j uf dr = 1, K = S, D and -E^ = k'^/M, and keeping only the first term, because the 
second term does not contribute to the Ml-transition. In this sense, our ^Si, ^Di-odd parity 
admixtures are only the part of the corresponding total wave functions which contribute to 
the PNC Ml-transition. 
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Figure|21(b) shows the odd parity ^Pq, ^Sq, ^Si and ^Z^i- admixtures for two potentials at 
AEry = 0.1 MeV. The ^Di-function is scaled additionally by \^Pd, where Pd is the probabil- 
ity of the D-state in a deuteron, because the corresponding Ml-transition is suppressed by 
this factor (P^*^"*^ = 0.0577, P^'^ = 0.0697). Again, one can see that in case of hard-core po- 
tentials, all wave functions vanish in the core-region, which leads to the relative suppression 
of the odd-parity ^Pq and ^So-components, whereas the ^Si and '^Di-configurations defined 
mainly by the long-range ttNN interaction are not sensitive to the potential at r > Tcore- In 
Fig. 121(c), we show the continuum wave functions at AE^ = 1 MeV. The main difference as 
compared with the previous case appears in the ^5*0 odd-parity admixture. It oscillates with 
r more strongly and has a node at r ~ 3.5 fm at AE^ = 1 MeV. This oscillating behaviour 
is manifested in the corresponding Ml-transition. 



IV. ASYMMETRIES 



The asymmetry of the deuteron disintegration in reaction with circularly polarized photon 
beam 



A 



RL 



o'x=i - cr\=-i 

(Ta=1 + (Ta=-1 



(22) 



consists of seven terms 



IT 



(23) 



i=l j=l 

defined by the interplay of dipole transitions caused by the parity-conserved and parity 
non-conserved interaction as follows 



^2 



T/7 
^3 



V2 



IT 



2 Re 

2 Re 

2 Re 

2 Re 

2 Re 

2 Re 

2 Re 
1 

2iV2 



T\Ml : D ^So) T{E1 : D '^Pq) /AT, 



T*(M1 : D ^So)T{El : 
T*{E1 : D 



'So) /AT, 



T*{E1 : D 
T*{E1 : D 
T*{E1 : D 
T*{E1 : D 
Tr [TT*]. 



^Po) T{M1 : D - 
3Pj) T(M1 : iPi 



3p 



^Pj) T(M1 : 3P] 



3p 



/A/", 
/A/", 



^Pi) r(Ml : D ^ 3Si) /A/", 
^Pi) T(Af 1 : D Hh)] /M. 



(24a) 
(24b) 
(24c) 
(24d) 
(24e) 
(24f) 
(24g) 
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Their explicit form in terms of the radial integrals read 



Re 



^2 



^3 



V2 



7 



7 



Re 



2 

3^3 AAM 

2 1 ji^ 



drru*{^Po : pr) [^(r) — A/2w(r)] 



• J drru*CSo : pr) m(^Pi : r) 
- V2/r) • J dru*CSo : pr) u{r) 
- • J dru*CPj : pr) ^I(^Pi : r) 



+ ) • / dru*CPj : pr) u^P, : r) 



7-2 * 



dru*{^Si : pr) ^(r) 



2 1 -3/2 
M 



Re 



1-2 * 



dru*{ Di : pr) i(;(r) 



(25a) 
(25b) 
(25c) 
(25d) 
(25e) 
(25f) 
(25g) 



Another asymmetry is related to the deuteron disintegration with unpolarized photon beam 
and polarized deuteron target 

An = ^^"=^ ~ ^^"=-^ , (26) 

CrMo=l + <^Mo=-l 

where M^, = 1(— 1) corresponds to the deuteron spin projection parallel (antiparallel) to the 
direction of the beam momentum. This asymmetry has also seven components 

= + (27) 



1=1 



Three of them, V^2,3j equal with the opposite sign to the corresponding — asymmetries 

= -V,\ = -V,\ Vi" = -VI (28) 



In these cases, the spin of the final states is zero and the corresponding Ml-transitions are 
proportional to 5_amo • The other four asymmetries are expressed as 



D 2 1 



Re 



- ■ j dru*CPj : pr) uCP, : 



r) 



2 1 
3jf 



Re 



M 



Us - 1/4 ^2* 



M 



dru*{^Pj : pr) uiy^Pi : r) 



--TT. 



2 ' 



TT. 



D 



2 ^ 



(29) 



The most important is the modification of 7rf. As we will see later, the spin-transitions in 
TTi and 7r2 proportional to a-re almost canceled in A^, but not m. A^- Therefore, the PNC- 
weak interaction of the vr-exchange may be clearly manifested only in the A^i-asymmetry. 
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V. RESULTS AND DISCUSSION 



We first discuss the A^-asymmetry. At — > i?thr, the Vi- and V2- terms only contribute 
to the total asymmetry. The signs of them are opposite and therefore their interference is 
destructive. The sign of the total asymmetry is defined by the dominant term. The strength 
of Vi^2 is determined by the values of the corresponding PNC-weak coupling constants and 
the behaviour of the proton-nucleon wave functions at short distances. For the case when the 
functions u{r) and u{^Sq : pr) are smooth at r < 1 fm (e.g. in the zero range approximation), 
one can neglect derivatives u' in Eqs. fll5|) and (fTTj) . Using the approximate expression for 
the Green function for r' < r and ~ 0: gi{E : r, r') ^ —r''^6{r — r')/3r, neglecting w and 
w', and taking into account the fact that the main contribution to the odd parity admixtures 
u{^Pi : r) and u{^Po : pr) comes from the terms proportional to f^{r'), one gets the following 
estimate 



V7 {hl-Jlhl){2 + x,) + hl{2 + x.) 
This estimation coincides with the result of the plane-wave Born approximation given in 



Ref. j2| and shows the dominance of the ^5*1 ^Pi PNC-transition with AJ = compared 
to the ^5*0 ^Pq with A/ = 0, 2. In case of the realistic NN-potential, the radial np-wave 
functions increase rapidly from zero at r = (for the hard core potential from r = r^ore) to 
the finite value at r ^ 1 fm. Since and \f'J\ decrease with r, the dominant contribution to 
the integrals in Eqs. (fTHj) and p7|l comes from the regions of r = 0.6 ~ 1.2 fm. This leads to 
increase of |V7/V2''| and to decrease of the asymmetries \Arl\ and |A_d|. Of course, we can 
not neglect the terms with derivatives u' because they are essential just in the region of the 
dominant contribution of the corresponding integrals. In our case u'{r), w'{r), u'{^So : r) at 
r < 1.2 fm are positive and large, especially for the hard-core (HJ)-potentials. In Eq. (fT3jl . 
the term proportional to u'{r) gives a constructive contribution and enhance IV2I, whereas 
in Eq. (|T7jl . u'^^ij/r) contributes destructively and suppresses |Vi|. As a result, we get the 
ratio of /V2 close to its raw estimate of Eq. (jHOI)- 

Figure a) shows the asymmetries A^l as a function of AE^ together with the partial 
asymmetries Vi and vTj. When AE^ increases, the PNC Ml transitions become important. 
At low Ai?^, asymmetries caused by the A/ = 0,2 PNC-forces and , generated 
by A/ = forces are close to each other numerically with the same sign. However, at 
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AE^/ ~ 0.5 MeV, decreases, changes sign and then its absolute value becomes much 
smaller than \ V2\, and it does not affect the asymmetry. In the limit of AE^ —>■ our result 
i^LR = 3.35 X 10^^) is in agreement with the previous calculations of the circular photon 
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polarization in the np — > reaction (P^ = (1.8 ~ 5.6) x 10" 

The PNC transitions with AI = 1 {AS = 0) are described by the irj, vTg and nj- 
terms, where terms are dominant and they are mostly determined by the weak vr-meson 
exchange interaction. In Fig. Efa), we show the vr^- asymmetry, the sum of ttJ + vTg -terms, 
and the coherent sum of all the AI = 1 transitions denoted as ir'^. At AE ~ 10 MeV, the 
absolute values of ttJ and vrj are the biggest among the other (Vi) terms and close to each 
other. But their signs are opposite. Therefore, the coherent sum is rather small 

<2 = ^J + ^2 - ^s{li - /£) ~ /X.O(Pd), (31) 



where and I\j are the radial integrals for the Ml-transitions in Eqs. ()25e|l and (j2| 
respectively. The finite value of tt12 is mainly caused by the non-symmetrical contribution of 
the deuteron d-wave in vr^, and t^I and it almost vanishes when = 0. In case of the zero 
range approximation in the limit AE^^ 0, this cancellation is exact In the real case 
the total contribution of the AI = 1 PNC interaction (vr''') is finite. However, its absolute 
value is smaller by a factor of 27 as compared with the result of Ref. 0]. Therefore, it 
seems to be difficult to get information about the AI = 1 PNC forces from 

The coherent interference of the V^- V2- and V^^-terms leads to sharp decrease of Aj^l 
down to zero at AE^ ~ 1.3 MeV (in case of Paris potential), and change a sign from 
positive to negative. Figure El (b) shows the total asymmetry Afi^ for the two potentials. 
For illustration, we also show the prediction of Ref. for the modified ZRA-model. One 
can see that the behaviour of the asymmetry A^l is similar qualitatively for the quite 
different models. In case of the HJ-potential, the asymmetry is smaller. The difference 
between two potentials at small AE^ = 0.01 ~ 1 MeV amounts to a factor of 2.5 ~ 3. The 
interce pt Af/r, = is shifted towards lower energies. The prediction of the modified ZRA- 
model 3] is close qualitatively to those of the Paris potential but the absolute value of A^l 
is much greater and the position of the intercept is shifted towards higher energies. This 
comparison with HJ-potential and ZRA-model has a rather illustrative character because the 
realistic potentials with the soft core repulsion are commonly accepted to be more adequate 
for description of the short range phenomena. From this point of view, only the prediction 
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obtained with Paris potential seems to be realistic. 

Figure |3fa) shows the A^j-asymmetry as a function of ^E^. There are two main differ- 
ences compared to the A^l asymmetry. Firstly, the components V2 and V4 are of the same 
sign. Secondly, there is no cancellation between the ^Pi ^Pj and D '^Si - transitions. 
Their coherent sum now behaves as 

vrf^ = vrf + vrf ~ - (32) 

and becomes a significant part of the asymmetry at large Ai?^. The sum of all transitions 
generated by the A/ = 1 PNC forces vr^ = 7r{^ + 7r|^ has the same sign as the V2- and 
V4- components. This leads to a non-monotonical behaviour of \A£)\ with a local minimum 
at Ai?^ ~ 2 MeV, but the sign oi A£, remains to be the same at < Ai?^ < 10 MeV and 
negative. In Fig.^J^b), we compare the results for Ad calculated with the two potentials. The 
difference between two asymmetries decreases with increasing the photon energy. However, 
the two results are similar in shape. 

The weak tt- meson exchange is mostly important at large ^E^. For illustration. Fig. 
shows the asymmetry Ad calculated as a function of /\E^ at different values of /it^ which 
cover its theoretical uncertainty: < /i^ < 2.5/1^*^*^*, where h}^^^^ is the "best" value of 
DDH. One can see that the constructive interference between week vr, and vector meson 
exchange results in increasing the absolute value oi Ad with increasing /i^r and leads to shift 
the position of the local minimum towards the lower energies. The absolute value of \Ad\ 
increases by a factor of 3 when R^^ changes from to 2.5 at 1 < AE^ < 10 MeV. 

Using the energy dependence of Arl and Ad-, one can obtain relations between the weak 
coupling constants. Thus, the standard representation of asymmetries through and /i^ 
read 

Arl = a^gphl + algphl + a° 5(^/1° + al{g^hl - gph^ + a'lgph'p^ + a^g^K, (33) 
Ad = hlgphl + h^gphl + hlg^hl + hl{gX^ - gph)) + h'^gph'^^ + Kg^K. (34) 

In the ideal case, having the asymmetries at six energy points and using the energy depen- 
dence of and hi one extract hi unambiguously. In practice, the number of "independent" 
equations for determination of hi is smaller, because some of {hi) are rather weak. The en- 
ergy dependence of the coefficients Oj and hi is shown in the Figs. IHl^a) and (b), respectively. 
For simplicity, we display only the dominant terms. 
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There are several points, where Arl and are particularly interesting. At AE — > 0, 
where the absolute values of both the asymmetries have a maximum, we get the following 
relations 

AnL ^ -(4.82^p/i° + 7A3gphl - 0.99^7^/i°) x 10-^ (35) 
Ad ^ -Arl. (36) 

. The point AE ~ 10 MeV can be used for analyzing the vr-meson exchange contribution in 
Ad: 

Ad - {lAQgphl - OMg^hl + 0.27gX^ - 0A3g^K) x 10^1 (37) 

The coefficient bj^ is governed by the long range interactions and therefore is not sensitive 
to the model of NN-interaction at short distances. 

The position of intercept Arl = at AE^ ~ 1.3 MeV may be also used for fixing the 
relation between coupling constants, but the experiment to find this position would be very 
difficult. On the other hand, another relations may be obtained when one of the term in 
Eqs. (jH^ and (jMj) vanishes but asymmetries have a finite and reasonable value. Thus, we 
have at AE^ ~ 0.4, a° = 0, and therefore 

Arl{^E^ ^ 0.4 MeV) ~ -{S.Ug.hl - O.Q7g^hl) x lO'l (38) 

Relations (jH^]) - (jHHj) are derived using the energy dependence of the coefficients a, and bt 
in Eqs. ()33|1 and (j34|l shown in Fig. 1^1 The later is defined by the short range behaviour of 
NN-forces, and is obtained with the Paris potential which has been, in particularly, designed 
for the adequate description various phenomena sensitive to the nucleon interaction at short 
distances. On the other hand, the Paris potential can not describe the neutron-proton 
scattering length which is its obvious disadvantage. Nevertheless, we convince that our 
results for the Paris potential would be coincide within ~ 20 — 30% accuracy with the 
predictions obtained with the other soft-core realistic potentials. This level of accuracy 
corresponds to the difference between our result and previous calculations of P^=Arl{E^ = 
Ethr) with different realistic potentials 
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VI. SUMMARY 



We have analyzed the energy dependence of two PNC-asymmetries in the deuteron photo- 
disintegration: one with circularly polarized photon beam [Aju) and another with polarized 
deuteron target (Ad). We show that by combining the measurements of A^xand A^, valu- 
able information on the PNC-nuclear forces may be obtained. Namely, using the energy 
dependence of Afi^ and A^, three constraints (equations) for determination of the PNC- 
coupling constants. 

Finally, we stress that the present investigation is a very first step. It would be impor- 
tant to verify if the predicted asymmetries are universal in the framework of other realistic 

The role of the 



potentials invoking the meson-exchange currents and relativistic effects js 
higher multipole transitions at higher energy is not quite clear. 

After completing this paper, the work by Liu, Hyun, and Desplanques has appeared 
in arXiv 3^. The authors have analyzed the A^j^- asymmetry using the realistic Argonne 
AVI 8 - potential. In spite of some difference in our models, the results of both papers are 
consistent to each other. Ref. 3^ gives A^ilIAE^ ~ 0) ~ +2.53 x 10~^, and Apu changes its 
sign at AE^ ~ 1.5 MeV. The contribution of the week vr-exchange transition is suppressed 
dynamically and it is about a factor of 30 smaller than the prediction of Ref. 2^. 



Acknowledgments 

We thank S. Date', H. Ejiri, C.-P. Liu, I. Khriplovich, and Y. Ohashi for fruitful discus- 
sion. One of author (A.I.T.) thanks M. Yasuoko, the director of Advanced Science Research 
Center, for his hospitality to stay at SPring-8. This work was supported in part by the 
Japan Society for the Promotion of Science (JSPS), and was strongly stimulated by a new 
project to produce a high-intensity MeV 7-rays by inverse Compton scattering at SPring-8. 



[1] E. M. Henley, Ann. Rev. Nucl. Sci. 19, 367 (1969). 

[2] E.G. Adelberger and W.C. Haxton, Ann. Rev. Nucl. Part. Sci. 35, 501 (1985); W. Haeberli 
and B. R. Holstein, in Symmetries and Fundamental Interactions in Nuclei, edited by W. C. 
Haxton and E. M. Henley (World Scientific, Singapore, 1996) pp. 17. 



17 



[3] B. Desplanques, J.F. Donoghue and B.R. Holstein, Ann. Phys. (N.Y.) 124, 449 (1980). 
[4] Ulf-G. Meifiner and H. Weigel, Phys. Lett. B 447, 1 (1999). 

[5] E. M. Henley, W.-Y.P. Hwang and L.S. Kisslinger, Phys. Lett. B 367, 21 (1996), and Phys. 

Lett. B 440, 449 (1998). 
[6] V.M. Dubovik and S.V. Zenkin, Ann. Phys. (N.Y.) 172, 100 (1986). 

[7] G. B. Feldman, G. A. Crawford, J. Dubach, and B. R. Holstein, Phys. Rev. C 43, 863 (1991). 
[8] W.C. Haxton, C.-P. Liu, and M. J. Ramsey-Musolf, Phys. Rev. C 65, 045502 (2002). 
[9] E.G. Adelberger, M.M. Hindi, CD. Hoyle, H.E. Swanson, R.D. Von Lintig, and W.C. Haxton, 
Phys. Rev. C 27, 2833 (1983). 
[10] S. A. Page et al., Phys. Rev. C 35, 1119 (1987). 

[11] C. S. Wood, S. C. Bennett, D. Cho, B. P. Masterson, J. L. Roberts, C. E. Tanner, and C. E. 

Wieman, Science 275, 1759 (1997). 
[12] V. V. Flambaum and D. W. Murray, Phys. Rev. C 56, 1641 (1997). 

[13] V.A. Knya'zkov, E.A. Kolomcnskii, V.M. Lobashov, V.A. Nazarenko, A.N. Pirozhkov, A.I. Sh- 
ablii, E.V. Shulgina, Yu.V. Sobolev, and A.L Yegorov, Nucl. Phys. A417, 209 (1984). 

[14] G.S. Danilov, Phys. Lett. B35, 579 (1971); G.S. Danilov, Sov. Nucl. Phys. 14, 443 (1972). 

[15] K.R. Lessey and B.H.J. McKellar, Phys. Rev. C 11, 349 (1975) and 12, 721 (E) (1975); M. 
Gari and J. Schlittcr, Phys. Lett. B59, 118 (1975); J.P. Leroy, J. Micheh, and D. Pignon, 
Nucl. Phys. A280, 377 (1977). 

[16] B. Desplanques, Nucl. Phys. A335, 147 (1980). 

[17] J.F. Cavaignac, B. Vignon, R. Wilson, Phys. Lett. B67, 148 (1977); 

[18] B. Desplanques, Phys. Lett. B 512, 305 (2001). 

[19] M. Snow et al., Nucl. Instrum. Methods, 440, 729 (2000). 

[20] W.-Y.P. Hwang, E.M. Henley, Ann. Phys. (N.Y.) 129, 47 (1980). 

[21] W.-Y.P. Hwang, E.M. Henley, and G.A.Miller, Ann. Phys. (N.Y.) 137, 378 (1981). 

[22] R. Schiavilla, J. Carlson, and M. Paris, Phys. Rev. C, 67, 032501 (2003). 

[23] B. Desplanques, Phys. Rep. 297, 1 (1998). 

[24] H.C. Lee, Phys. Rev. Lett. 41, 843 (1978). 

[25] T. Oka, Phys. Rev. D 27, 523 (1983) 

[26] LB. Khriplovich and R.V. Korkin, Nucl. Phys. A690, 610 (2001). 
[27] U. G. Meissner, Mod. Phys. Lett. A 5, 1703 (1990). 

18 



n 



TABLE I: Weak coupling constants determined from the "best value" of Ref. j^]. All values are 
given in units of 10~^ 





h\ 


" p 


h? 
"'p 




hi 


K 


-1.14 


-0.02 


-0.07 


-0.95 


-0.19 


-0.11 


0.46 



[28] R.V. Korkin, arXiv: nucl-tli/0206044| 

[29] M. Lacombe, B. Loiseau, R. Vinh Mau, J. Cote, P. Pires, R. de Tourreil, Phys. Lett. B 101, 
139 (1981). 

[30] M. Lacombe, B. Loiseau, J.M. Richard, R. Vinh Mau, J. Cote, P. Pires and R. de Tourreil, 

Phys. Rev. C 21, 861 (1980). 
[31] T. Hamada and I.D. Johnston, Nucl. Phys. 34, 382 (1962). 

[32] D.B. Kaplan, M.J. Savage, R.P. Springer, M.B. Wise, Phys. Lett. B449, 1 (1999). 
[33] M.J. Savage, R.P. Springer, Nucl. Phys. A686, 413 (2001). 

[34] C.A. Barnes, J.H. Carver, G.H. Stafford, and D.H. Wilkinson, Phys. Rev. 86, 359 (1952). 

[35] Y. Birenbaum, S. Kahane, and R. Moreh, Phys. Rev. C 32, 1825 (1985). 

[36] W. Tornow, N.G. Czakon, C.R. Howell, A. Hutcheson, J.H. Kelley, V.N. Litvinenko, S.F. 

Mikhailov, I.V. Pinayev, G.J. Weisel, H. Witala, Phys. Lett. B574, 8 (2003). 
[37] E.M. Henley, Nucl. Phys. A300, 273 (1978). 

[38] J. Carbonell, B. Desplanques, V.A. Karmanov, and J.F. Mathiot, Phys. Rep. 300, 215 (1998). 
[39] C.-P. Liu, C.H. Hyun, and B. Desplanques, |micT-th/0403009, 



19 



E 




E 



10' 



10^ 



10' 



10' 



yD^np 






Paris 




1-1 1 ■ — - 

HJ 




(b) . 



2 4 6 8 10 
E - E„, (MeV) 



2 4 6 8 10 
E - E„, (MeV) 



FIG. 1: The total cross section of the deuteron photo-disintegration as a function of the energy 
excess AE^ = — -Ethr- (a) Result for the Paris potential. Contributions of the Ml and E\ 
transitions are shown by the dashed and dot-dashed curves, respectively, (b) The total cross 
section for the Paris (solid) and Hamada-Johnston (dashed) potentials. The experimental data on 
the total cross section are taken from Refs. 3^ (open 



Ml-transition (filled squares) are taken from Ref. 
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FIG. 2: The odd-parity admixture to the proton-neutron wave functions calculated with the Paris 
(solid curves) and HJ (dashed curves) potentials, (a) Results for the deuteron wave functions, (b) 
and (c) Results for the continuum np- wave functions at /S.E^ = 0.1 and 1 MeV, respectively. 
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ErE,.r (MeV) E -E„, (MeV) 

FIG. 3: Asymmetry of the deuteron disintegration in the reaction 7D pn with circular polarized 
photon and unpolarized deuteron as a function of energy excess E-y — Ethr- (a) Relative contribution 
of the different odd-parity transitions for the Paris potential. The sign in the bracket denotes the 
sign of the corresponding term, (b) Comparison of the total asymmetry for the Paris (solid), 
Hamada- Johnston (dashed) potentials and the modified ZRA of Ref . j22i] (dot-dashed) . 




ErE.,r (MeV) E -E,,, (MeV) 

FIG. 4: Asymmetry of the deuteron disintegration in the 7-D — > pn reaction with polarized deuteron 
and unpolarized photon beam as a function of energy excess -E^ — -Ethr- (a) Relative contribution 
of different odd-parity transitions for the Paris potential. Notation is the same as in Fig. |31 (a), 
(b) Comparison of the asymmetry for the Paris and Hamada-Johnston potentials. 
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FIG. 5: Asymmetry of the deuteron disintegration in the — > pn reaction (j4d) with different 
values of the PNC 7r-exchange coupUng constant: R = fn/fn^^^ = 0, 1, 2.5, where f^'^^^ is the "best 
value" of Ref. 




FIG. 6: (a) The quantities of Eq. (|35|) . (b) The quantities bi of Eq. (|36|) . We display only the 
largest components. Results are obtained with the Paris potential. 
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